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The bound state and low-energy scattering properties of two oriented dipoles are investigated for 
both bosonic and fermionic symmetries. Interestingly, a universal scaling emerges for the expectation 
value of the angular momentum for deeply-bound two-dipole states. This scaling traces to the 
pendulum motion of two dipoles in strong dipole regime. The low-energy scattering phase shifts 
of two dipoles also show universal behavior. These universal observations make connections to the 
scaling laws reported in Refs. [1, 2] for three dipoles. Atomic units are used throughout this work. 

PACS numbers: 



I. INTRODUCTION 

The recent experimental realization of ultracold po- 
lar molecules [3] has attracted tremendous attention 
through its potential future applications in studies of 
astrophysics, condensed matter physics, quantum com- 
puting, and ultracold chemistry [4, 5]. It also pro- 
vides great opportunity for studies of few-body physics 
due to the long-range, anisotropic nature of the molec- 
ular interactions in the presence of an external electric 
field. For polar molecules that have exothermal reaction 
paths [6] , the high short-range reaction probability leads 
to large ultracold reaction rates that have been observed 
in experiments and have also been explained by a uni- 
versal theory [7-13]. Novel phenomena involving polar 
molecules [4, 5] are in particular expected from the long- 
range interactions. When the molecules are nonreactive, 
rich resonant features that are tunable via external field 
have been predicted in the scattering of two molecules 
in different geometries [14-23]. Interestingly, these res- 
onances are formed with different mechanisms [17] and 
therefore occur in an irregular order. 

Whereas two-dipole physics has proven to be surpris- 
ingly complicated, the threc-dipole physics near two- 
dipole resonances manifests simple universal behavior for 
both bosonic [1] and fermionic [2] dipoles. This is quite 
counterintuitive considering that a three-body problem 
is generally more complicated when compared to a two- 
body problem. Moreover, in contrast to the previously 
known three-body systems where non-universal short- 
range physics is important even at unitarity [25] , a system 
with three dipoles at unitarity can be universally defined 
by two-dipole physics that is characterized by the dipole 
length defined as 



de = md^/2, 



(1) 



where dm is the magnitude of the electric dipole mo- 
ment, and m is the mass of a dipole. A close relation- 
ship between two- and three-dipole problem is also man- 
ifested in the effective repulsion between a dipole and a 
deeply-bound dipolar dimer that has been identified in 
Refs. [1, 2]. 

In the present work we study universal two-dipole 
physics that can be relevant to the three-dipole physics 



mentioned above. In particular, the properties of dceply- 
and weakly-bound dipolar states as well as the behavior 
of low-energy scattering phase shift are discussed for both 
bosonic and fermionic dipoles. It will be shown that when 
the dipolar interaction is off-resonant, the binding energy 
of two dipoles scales like 1/md'j, and the size of the corre- 
sponding state grows linearly with df . For deeply-bound 
states, the expectation value (L^) shows universal \/di 
scaling, where L is the angular momentum operator. As 
will be shown, this scaling traces to the pendulum motion 
of two dipoles in a strong external field. Finally, universal 
scaling behavior will be shown for both the real part and 
and the imaginary part of the phase shifts that charac- 
terize elastic scattering and the scattering into different 
angular momentum channels, respectively. 

This paper is organized as follows. Section II intro- 
duces our model and method for solving the two-dipole 
problem. Section III discusses the bound state proper- 
ties of the dipoles, and Section IV discusses the scat- 
tering properties. Finally, we summarize our results in 
Section V. 



II. THEORY 

To a good approximation at long range the dipolar 
molecules in an external field can be treated as point 
dipoles fully-oriented along the field direction. The 
Schrodinger equation for two dipoles in spherical coor- 
dinates is given by 



m r dr"^ 



(2) 



where Vdd is the interaction between two dipoles that are 
oriented to the direction of the external field z: 



dd 



2di 1 - 3(z • f)^ 



fir). 



(3) 



The function /(r) = tanh(r/ro)^^ cuts off the dipolar 
interaction around the short-range length scale tq, which 
avoids an unphysical collapse of the system in the case of 
the — 1 /r^ singularity at the origin. 
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A. Diabatic representation 

To solve Eq. (2), it is natural to expand ip in the basis 
of angular momentum cigcnstatcs \lmi) as 



(4) 



where different to; decouple when the quantization axis is 
along the direction of the external field. For large r, dif- 
ferent \lmi) states decouple and physically serve as scat- 
tering channels for two dipoles. In this basis, Eq. (2) 
reduces to a set of coupled equations: 



1 + 
TO dr'^ mr^ 



Fp{r)+ 



de. 



(5) 



where the coupling coefficients D^{mi] I, I') are expressed 
in terms of 3-j symbols: 



D^imrJJ') ==(_ip+i4V(2? + l)(2?' + l) 
I 2 l' \ f I 2 1' 







-TO/ TO; 



(6) 



Note that I and V are coupled only when |/' — ^| ^ 2, and 
specifically when |/' — Z| = 2 for indistinguishable dipoles 
that are considered in the present work. 



B. Adiabatic representation 

Another way to solve Eq. (2) is to use the adiabatic 
representation, as has been implemented in Refs. [VI- 
20]. It has been shown [17-20] that the adiabatic rep- 
resentation gives better characterization of the various 
dipolar resonances [17, 19, 20], presumably because the 
adiabatic channels are less coupled than the angular mo- 
mentum basis at small distances. At large distances the 
adiabatic basis reduces to the angular momentum basis, 
so both representations will give the same scattering ma- 
trix. Although these two representations are physically 
equivalent, we would like to compare their efficiencies and 
use the one that is the most convenient for our study that 
involves both analytical and numerical work. 

In the adiabatic representation, the two-dipole wave- 
function is expanded in the adiabatic basis (pi, as 



r 



(7) 



where Q, represents the polar angle 9 and the azimuthal 
angle Lp, v and v' are the adiabatic channel indices. The 
adiabatic channel functions 0™' are obtained by solving 
the adiabatic equation 

( ^ + 14d 1 C (r ^) - f/.(r)C (8) 



Upon substitution of the adiabatic expansion Eq. (7), 
Eq. (2) also takes a multi-channel form: 



TO dr^ 



Frir) + 



771 ^ ^ 



dr 



dr 



-P. 



F.™' (r) 



-Y.Q.^..Fl^^{r)^EF:-^^{r), (9) 



but is now characterized by the non-adiabatic couplings 
Pu,v'{r) and Qv.v'i^)^ defined as 



d 
dr 



— f 

dr ' 



(10) 

The double brackets in the above definition indicate in- 
tegration over the angles ft. 

In the asymptotic region (r 3> di), the adiabatic chan- 
nel functions 4>i, can be calculated perturbatively from 
\lmi), which allows us to analytically derive the asymp- 
totic form of the non-adiabatic couplings. To calculate 
the leading order of the non-adiabatic couplings Pi,y for 
channels v and v' that asymptotically approach the di- 
abatic channels {l,^,mi) and {I,,', mi), respectively, the 
following expansion of the channels functions is required: 



(11) 



The expansion coefficients ry^ are obtained by \n — vl-th 
order perturbation theory: 



where 



V:^C:Adi/rr-'\ (12) 
D3{mi;l, + 2kJ, + 2k±2) 



n—v^l 

n i;(/, + l)-(/, + 2fc±2)(L + 2/c±2 + l)' 



C^, =1. 



(13) 



In the above expression the upper sign is taken when 
n > V, and the lower sign, n < v. The leading order term 
in the asymptotic expansion of P^.^' is then calculated as 



p , 



i\v-v'\ v'^l 



n—u 



\v'-n\ClCi. (14) 



Here the upper and lower signs are taken for v' > v and 
v' < v, respectively. 

Deriving the leading order expression for Qj^y requires 
an expansion of (p^, with two more terms: 



(15) 



where the upper and lower signs are taken for v' > u and 
v' < V, respectively. This expansion gives, for v' — 



g,,.^^[(C)2 + (c^V)' 



(16) 
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FIG. 1: (color online) Number of channels required for con- 
verging the dipolar bound states within 1% in the diabatic 
(spherical harmonic) and the adiabatic representations. 



for — iy\ = 1 



y^v^x) ^v'±\^ y^v'^\) ^v^\ 



}; (17) 



and for \v' — i^l ^ 2, 



.\l^-l^'\+4: 



n—uzkl 



(18) 



In all the above expressions for Qu,v' , the upper and lower 
signs are taken for v' > v and i'' < v, respectively. All 
the above asymptotic behavior for the non-adiabatic cou- 
plings has been verified by our numerical calculations. 



C. Comparison of diabatic and adiabatic 
representations 

In the following we give a brief comparison of the dia- 
batic and adiabatic representations concerning their con- 
vergence with number of channels. Figure 1 shows the 
convergence pattern for the energies of the ground state 
and the highest excited state with respect to the number 
of channels as increases. Generally, the adiabatic rep- 
resentation converges more quickly, and particularly for 
deeply-lying states it gives better convergence for larger 
dg. For highly- lying states more channels are required in 
order to reach convergence for larger dg in both represen- 
tations. 

The quick convergence of the adiabatic representa- 
tion for the deeply-bound dipolar states makes it most 
advantageous in our analytical study of these states in 
Section III A. For studies of the weakly-bound states or 
scattering, however, the adiabatic potentials and non- 
adiabatic couplings for high-lying channels at small dis- 
tances (r < df) cannot be analytically determined and re- 
quire numerical diagonalization. Therefore even though 



the adiabatic representation gives faster convergence it 
does not provide significantly better performance, and we 
will use the diabatic representation since it can be more 
easily generated in calculations of weakly-bound states 
and scattering. 



III. BOUND STATE PROPERTIES 
A. Deeply-bound states 

In the presence of a strong external field, dipoles in a 
deeply-bound state tend to align themselves in the lin- 
ear configuration where the dipolar potential Eq. (3) is 
angularly minimal. For small polar angles 9, Eq. (3) can 
be approximated by 



2d, 



302 



(19) 



which corresponds to a harmonic potential in the 6 di- 
rection with trapping frequency of a; = y^2Adi/m'^r, and 
a zero-point angle given by: 



muj \ r 



-1/4 



(20) 



As di increases, two dipoles in a deeply-bound state 
become more angularly localized and undergo pendu- 
lum motion within an angle roughly determined by 6q. 
The expectation value of the angular momentum (L^) is 
therefore expected to grow with d,. 

In view of the connection between (L^) and the barrier 
in the collision of a dipole and a dipolar dimcr [1, 2], the 
following studies the scaling of {L^) with dg. The adia- 
batic approximation of the two-dipole solution, Eq. (7), 
takes the following form: 



mi 
n.u 



(21) 



where n indicates the approximate quantum number for 
the radial motion. The angular wavefunction 0™' for the 
harmonic potential Eq. (19) can be directly written down 
in terms of the associated Laguerre polynomials L"(x): 



2 r(;^ + l) 
elV{u+\mi\ + l) \0o 

1 



"H I / /)2 



2 flT 

X e '^o 



27r 



(22) 



By using the approximate angular momentum operator 
for 6* < 1: 



i2 



I d B 1 ^2 



'86 09 e'^difi^ 



(L^) can be determined as: 



(23) 



(24) 
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FIG. 2: (color online) The scaled angular momentum ex- 
pectation value {L^) /{di/roY^'^ for bosonic dipoles (a) and 
fermionic dipoles (b). 



Since for deeply-bound states the radial wavefuction F™^ 
is localized around r <r[)^ the radial integral should scale 



with 1/ 
scaling 



ro and (L ) is expected to obey the following 



(25) 



It is interesting to note that this scaling is independent 
of the identical particle symmetry. This can be simply 
understood from the angular localization of a pendulum 
wavefunction that makes the exchange effect negligible. 

Figure 2 shows the angular momentum expectation 
value (L^) calculated numerically. The quick changes in 
(L^) at small dg is related to the transition of the dipolar 
states from a single angular momentum character to a 
pendulum character. The angular and radial excitations 
for the corresponding dipolar states can be seen from 
the wavefunctions shown in Fig. 3. The equally spaced 
(i^) from numerical calculations suggests the validity of 
the approximated scaling in Eq. (25). The radial excita- 
tions n introduce splittings on top of the scaling given by 
Eq. (25), but such splittings become finer with a sharper 
short-range cutoff function /(r) as both expected from 
Eq. (24) and verified by our numerical calculations. 



B. Weakly-bound states 




FIG. 3: (color online) The cuts of the deeply-bound two- 
dipole wavefunction along the x-z plane for bosons. Differ- 
ent angular and radial excitations {\rai\,v,n) are shown: (a) 
(0,0,0), (b) (1,0,0), (c) (0,1,0), (d) (0,0,1), (e) (1,0,1), (f) 
(0, 1, 1). The dipole length di = 600ro for all these wavefunc- 
tions. 
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FIG. 4: (color online) The characteristic size Vc (a) and 
the characteristic energy \Ec\ (see text) (b) for both bosonic 
and fermionic dipoles with |m; |=0 and 1. The irregular jumps 
come from the avoided crossings between states with different 
angular momentum character. 



For weakly-bound states, we first discuss the proper- 
ties of such states away from a dipole-dipole resonance. 
These properties are important in determining the scal- 
ing laws for three-dipole recombination [1] and possi- 
bly also for the vibrational relaxation of a weakly-bound 
dipolar dimer in an off-resonant situation. To this end, 
we introduce the characteristic size rc and the character- 
istic energy Ec for the weakly-bound dipolar states, as 
the expectation value of r and the energy respectively for 
the state that is right below a zero-energy bound state. 



Figure 4 shows numerically calculated rc and Ec as 
de increases. In all cases it is shown that Vc oc de and 
Ec oc 1/md'j. These scalings contrast interestingly with 
the deeply-bound states, where the sizes of the states 
shrink to rg and the energies grow deeper linearly with 
di- 

Near a dipole-dipole resonance where a zero-energy 
bound state is formed, the binding energy of a weakly- 
bound state is closely related to the low-energy expansion 
of the scattering phase shift. In the presence of s-wave 
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10"^ 
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single-channel scattering with a potential [26-28]: 
Re[(5™'(fc)] ^ —a^'k, where fc = \JmE is the scattering 
wavenumber. 

The effective range expansion for short-range poten- 
tials [29, 30] has the following low-energy expansion of 
the phase shift: 



(27) 



where the power of fc increases by 2 for consecutive 
terms. For dipolar scattering, however, our numerical 
study show that the phase shifts are in expansions with 
increment of fc for all partial waves: 



FIG. 5: (color online) The scaling behavior of the bind- 
ing energy Eh with di and V;"'' near p-wave resonances for 
fermionic dipoles. The dipole length di is tuned to have one 
or two bound states (BS) for m; = and mi = 1. 



contribution (mi—0 for bosons) the binding energy Eb 
can be written in terms of the s-wave scattering length 
as El, « l/"T-(a™o)^ finding the pole in the scat- 
tering amplitude [24]. When s-wave contribution is ab- 
sent, searching for a pole in the scattering amplitude is 
not straightforward due to an anomalous low-energy ex- 
pansion of the phase shifts [see Eq. (30)]. Nevertheless, 
near a p-wave resonance for fermions we have numeri- 
cally identified that Eb cx d^/V;™{, where the scattering 



Re[C(fc)] = -«rfc 



Vpk^ + 0{k^). (28) 



volume V^™' diverges near a p-wave resonance. As shown 
in Fig. 5, the proportionality constant in the scaling of Ef, 
depends on mi but is independent of the short-range in- 
teraction details as the number of bound states changes. 



Our discussion here will be restricted up to the term 
that starts showing non-universal behavior. As has been 
shown in previous works [14-20], the scattering length 
a™' shows non-universal resonant behavior for ^=0. We 
have further verified that in this case the resonant be- 
havior persists in all higher terms. For fermionic dipoles, 
our numerical study shows non-universal resonant behav- 
ior starting from V^"'|, while b]^^ remains universal. Non- 
universal resonant behavior is therefore expected to start 
from the term of fc^'^^^ for only I—Iq, where is the low- 
est partial wave allowed for a given symmetry. For the 
terms lower than fc^'""*"^ the coeflticients are expected to 
be universally determined by d£. 

For l>0, Ref. [18] has analytically derived a universal 
expression for the T-matrix to the leading order in k, 
their results can be readily used to give the following 
expression for al"' , 



IV. LOW-ENERGY SCATTERING 
PROPERTIES 

In previous studies of two-dipole elastic scattering [14- 
23], the behavior of the scattering cross-sections have 
been studied quite extensively. Our goal here, however, 
is to study the two-dipole scattering by examining the 
low-energy expansion of the phase shifts. Of particu- 
lar interest are the phase shifts from the lowest partial 
wave, which is the most sensitive to the dipole-dipole 
resonances at zero energy. 

In the present case where the angular momentum is 
not conserved, we define the phase shifts 5^' from the 
diagonal scattering matrix elements 



Ll 



HSZ')/2^. 



(26) 



In the presence of couplings between different partial 
waves, (5["' acquires an imaginary part that character- 
izes off-diagonal scattering amplitudes. The real part of 
controls the elastic scattering. While dipole-dipole 
scattering is multichannel, some insights into the low- 
energy expansion for Re[(5™'(fc)] can be gained from the 



D^{mi\l,l) 
' 2l{l + l) ■ 



(29) 



Next we study the scaling behavior of 6™' and the 
scattering volume V;™' by numerical calculations. The 
universality is tested by adding an isotropic short-range 
interaction Mso = Vbsech^(r/ro). Figure 6 (a) shows the 
di dependence of a few with different Vq. It is clearly 
seen that 6™' follows a scaling behavior that is inde- 
pendent of Vq. An exception for this universal behav- 
ior is h^Q^, where the non- universal resonant behavior 



already begins in the lower term a 



mi=0 
1=0 



Nevertheless, 



is found to follow a universal df background scal- 
ing with non-universal resonant features on top. Table I 
lists some numerically determined scaling coefficients. 

The study of scattering volume is more challenging due 
to the difficulty to numerically fit Eq. (30) to the third 
order accurately. Our numerical study shows that the 
scattering volume Vj'™'' follows a df scaling in general. 
For 1=1, the lowest partial wave allowed for fermionic 
scattering, non-universal resonant features are expected 

for Vi2i^- Nevertheless, a universal d| background scal- 
ing can be identified from Fig. 6 (b). The positions of 
the resonant features clearly depend on the short-range 
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mi 



-2.60x10"^ 
2.56x10"^ 
2.30x10"^ 
-8.32x10"^ 
-5.70x10"^ 
1.27x10"^ 



TABLE I: The numerically calculated universal scaling for the 
phase shift expansion coefficient bj"' for a few symmetries. 
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FIG. 6: (color online) (a) The d^-dependence of the coeffi- 
cient in the low-energy expansion of Re[5™' (fc)]. (b) The 
d^-dependence of the scattering volume V;™' . 



interaction as tuned by Vq, but the background scaling 
remains unaltered. 

Finally we discuss the low-energy expansion for the 
imaginary part for the phase shift Im[(5™'(fc)]. The fol- 
lowing expansion 



Im[5p(fc)] 



-cj"'fc2-|-0(fc3) 



(30) 



is found from our numerical calculations, indicating a 
vanishing imaginary part in the scattering length when 



kdi <C 1. By using the unitary constraint on the diago- 
nal element of the 5'-matrix to leading order, c™' can be 
determined as 



1 
12 



D:,{mi-l,l-2) 



for all partial waves. 



J3(m;;/,/ + 2) 
{l + l)(l + 2) 



(31) 



V. SUMMARY 

To summarize, we have studied the universal proper- 
ties for two dipoles. The long-range, anisotropic dipo- 
lar interaction brings rich, universal physics that has 
key implications for the universal three-dipole physics. 
This is shown particularly for both the deeply-bound 
and weakly-bound sides of the dipolar spectrum. For 
the deeply-bound dipolar states, the pendulum motion 
between the dipoles gives rise to a universal growth in 
the expectation value of the angular momentum, which 
produces a centrifugal barrier between a dipole and a 
dipolar dimer [1, 2]. For the weakly-bound states, gen- 
eral scalings of the binding energy and the size of the 
states are identified, despite the complicated level cross- 
ings for states with different angular momentum charac- 
ters. Finally, the low-energy scattering phase shifts for 
two dipoles are predominantly determined by the dipole 
length, with some non-universal ingredients that give rise 
to resonant features. 
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